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1. Introduction 

Form factors measured in electromagnetic and weak processes are fundamental probes of 
hadron structure. Calculations of such observables using lattice QCD and, in particular, the nu- 
cleon form factors [[I], ||, ||] has intensified during the last couple of years due to improvements 
which allow full lattice QCD calculations with controlled lattice systematics [Q]. The focus of the 
current work is the study of the electro-magnetic (EM) and weak N to A transition form factors 
(FFs). Experiments on the N to A EM transition have yielded accurate results on the EM transition 
form factor for low momentum transfer []|] that point to deformation of the N/A system. The axial 
N to A transition FFs are experimentally not well known but there are ongoing experiments using 
electroproduction of the A resonance to measure the parity violating asymmetry in N to A. Lattice 
QCD enables calculation of these fundamental quantities from first principle. Our previous calcu- 
lation of these form factors utilized quenched and dynamical Wilson as well as a hybrid scheme 
with domain wall (DWF) valence quarks on an improved staggered sea []6[ [7|, |8p. A study of the 
N to A transition using chiral dynamical quarks in a unitary approach is presented in this work 
where, in addition, we employ the coherent sink method [Q] in order to achieve the better statistical 
accuracy on the determination of the form factors. 

2. Lattice Techniques 

We use Nf = 2 + 1 dynamical domain wall fermions generated by the RBC and UKQCD 
collaborations [0]. The lattice spacing a -1 = 1 .73(3) GeV is fixed using the Q. mass. The length of 
the fifth dimension is taken sufficiently large to suppress chiral symmetry breaking. Fixing L^/a = 
16 gives an additive residual mass ~ 10% of the light quark mass used in this work. We consider 
configurations on a lattice of size 24 3 x 64 corresponding to pion mass of 0.331(1) GeV. We use the 
standard interpolating operators to create nucleon and A states and employ gauge invariant gaussian 
smearing of the quark fields with APE-smeared gauge fields optimized for best suppression of 
excited states for the nucleon [QJ. Suppressing excited state contributions in the three-point function 
is particularly crucial since for this study a source-sink separation of 0.9 fm is used. We show in 
Fig. |] that extending the source-sink separation to 1.14 fm the plateau values for the dominant 
dipole form factor Gmu which are the most accurate, are consistent with a time-separation of 
0.9 fm, but with a two-fold increase in statistical errors. 

The three-point functions that are needed are given by 

{G^\t 2 ,tu\>\\>;r T ))= £ e-^^e + ^r p t a (a\T[xr(^2WxiA)x p N M} |G> (2-1) 

X 2 ,Xl 

where (jc) is a local current, q = p' — p is the momentum transfer, a is the Lorentz vector index 
for the A and T T projection matrices in Dirac space [^J. The large Euclidean time limit of the ratio 
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yields a time-independent function n ff (p',p ;r T ;/i) (plateau region). In addition, all field renor- 
malization constants cancel and therefore n CT is a combination of the Lorentz invariant form factors 
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and known kinematical factors. We use sequential inversions through the sink to evaluate the three- 
point function of Eq. (|z~l|). In this method the quantum numbers of the hadron are fixed, which 
means that a particular value of a and r T must be chosen. This freedom is exploited in the construc- 
tion of sources for the sequential propagator with the goal to produce optimal linear combinations 
of n CT involving a maximal set of momentum vectors, thereby obtaining a maximum number of 
statistically independent measurements It turns out that three such sinks suffice for achieving 
this goal and enable us to extract the momentum dependence of the electromagnetic, axial and 
pseudoscalar N to A FFs accurately. A new ingredient of the current work is the use of the coherent 
sink technique [Q] in order to reduce the statistical noise. This consists of creating four sets of 
forward propagators for each configuration by placing sources at: 
(0,0), (L/2,16), (0,32) and (L/2,48). 

From each source (x,-,7}), a zero-momentum projected A source is constructed at To away, i.e. at 
(xj,Tj + T a ) and a single coherent backward propagator is calculated in the simultaneous presence 
of all four sources. The cross terms that arise vanish by gauge invariance when averaged over 
the ensemble. The forward propagators are already computed by the LHPC Collaboration [§] and 
therefore we effectively obtain four measurements at the cost of one. This assumes large enough 
time-separation between the four sources to suppress contamination among them. An open ques- 
tion is whether there exists correlation among these four measurements. In Fig. [2| we show the 
dependence of the jackknife error on Gm\ for different coherent sink bin sizes, which verifies that 
cross-correlations between the different sinks are absent. 

The full set of data obtained at a given Q 2 value is analyzed simultaneously by a global % 2 
minimization using the singular value decomposition of an overconstrained linear system All 
the results presented here are obtained by analyzing 200 configurations or a total of 200 x 4 = 800 



measurements of the ratio given in Eq. (2.2) 
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Figure 1: The ratio Si of Eq. (3.2) versus t/a for 

a source-sink separation 0.91 fm shifted by a time Figure 2: Dependence of the jackknife error for 
slice (blue triangles) and 1.14 fm (red circles) for the G M \ (Q 2 ) on the coherent sink bin sizes, 
smaller non-zero q 2 . 
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3. Electromagnetic N to A Transition form factors 

The electromagnetic transition matrix element is decomposed in terms of three Sachs (FFs) 

(A(p', S ')\j fl \N(p,s)} = i^ ( ^-^^ y 2 u a {p'y)ff atlU {p,s) (3.1) 

with 

6c» = G M i (<7 2 )<! + G E2 {q 2 )K E a l + G cl (q 2 )K% 

where K^,K E2 and are known kinematical factors In this work we present results for 
the dominant magnetic dipole form factor Gm\{<J 2 )- Following Ref. [Q] we construct the optimized 
three-point function S\ from which Gmi(Q 2 ) is directly determined 

5i(q;/i) = £ n a (Q,-q;r 4 ;ix)=iA^(p 2 -p3)8 Ufi + (p3-pi)52^ + (pi-p2)8 3 ^GMi(Q 2 ) ■ 

(3.2) 
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Figure 3: Gm\{Q 2 ) using DWF fermions and using 
the hybrid action. The diamonds show experimental 
results. The solid (dashed) line is a fit to dipole (ex- 
ponential) form for the DWF data. 
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Figure 4: Cj for DWF, the hybrid action and 
quenched Wilson fermions ( m n =410 MeV) B. 
The green line is a dipole fit to experimental 
data JlOtl ■ The solid (dotted) line is a fit to dipole 
(exponential) form of the DWF data. 



In Fig. U we show the results of this work on Gmi(Q 2 ) using DWF. These are compared 
with previous results obtained with a hybrid action that uses Asqtad improved staggered fermions 
generated by the MILC collaboration and domain wall valence quarks [0]. The pion mass in the 
DWF calculation is 33 1 MeV and in the hybrid action 350 MeV. These values are close enough to 
allow a direct comparison. Indeed the results are in very good agreement. Fits to a dipole form, 
go/(l + 2 2 / m o) 2 ' as we U as to an exponential form go exp(— Q 2 /riio 2 ) described equally well the 
lattice results. A compilation of the experimentally available data (for more details see Ref. [{7]]) is 
also shown in Fig. || showing a clear disagreement between lattice results and experiment. This is 
reflected in the value of the dipole mass of mo = 0.78 GeV obtained by performing a dipole form 
fit to the experimental data as compared to mo = 1.164(20) GeV for the lattice results. A possible 
explanation for the faster falloff of the experimental data maybe the lack of significant chiral quark 
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effects -or equally the lack of strong pion cloud- from the still heavy pion mass ensembles that are 
utilized. Similar behavior is also observed for the nucleon electromagnetic form factors [jl]], that 
may again point to the importance of chiral quark effects. The N to A case is particularly clean 
since there is no ambiguity regarding disconnected contributions and thus the flatter dependence 
observed in the N to A EM FFs must be of different origin. The large disagreement observed here, 
however, would require large pion cloud effects to set in as we lower the pion mass. Such large 
pion effects have been shown to arise in chiral expansions [|ll|] and it is thus interesting to perform 
the calculation for m K < 250 MeV where they are expected to set in. We are currently analyzing 
results to extract the subdominant FFs, Gei and Gqi using the same DWF configurations. 



4. Electroweak N to A Transition form factors and Goldberger-Treiman relations 

We consider nucleon to A matrix elements of the axial and pseudoscalar currents defined by 

x a % a 
A^(x) = \j/(x)Y^Y 5 — \ l /(x) , P a (x) = \j/(x)y 5 — y(x) (4.1) 

where r" are the three Pauli-matrices acting in flavor space and y the isospin doublet quark field. 
The invariant proton to A + weak matrix element is expressed in terms of four transition form factors 
in the Adler representation as 

<A(i/y)W^)> = *vl(^w) 1/4 *^) 

fC^(q 2 ) Ci(q 2 ) lv \ ( \P_ur*f^ 

7 + — P [gXiigpv-gXpgnv)q P +C 5 {q )g Xll + 2 qxqp 

V m N / m N 



u P (p,s). (4.2) 



The form factors C3 (q 2 ) and C^(q 2 ) belong to the transverse part of the axial current and are both 
suppressed [||] relative to the dominant form factors C^{q 2 ) and C^(q 2 ). The latter two are the 
equivalent to the nucleon axial FFs Ga(Q 2 ) and G P {Q 2 ) respectively [0]. 
The pseudoscalar transition form factor G n N&{q 2 ), is defined via 

2m q <A(p,s)\P pv(p,,)> =l ^ __ j — — u A+ (p,s)—u P {p,s) .(4.3) 

Taking matrix elements of the axial Ward-Takahashi identity d^A^ = 2m q P a leads to the non- 
diagonal Goldberger-Treiman (GT) relation 

C$tf) + £&tf) = ± G ™¥ )f f< . (4.4) 
m l N 2m N m l n - q 1 

The PCAC relation on the hadronic level d^A"^ = f^m^Ti", relates the pseudoscalar current to the 
pion field operator and therefore provides the connection to the phenomenological TiNA strong 
coupling g K NA = GWa(O) that appears in Eq. (AA). Assuming pion pole dominance we can relate 
the form factor Cg to G^na via: 

1 r^f^\ 1 G nNA {q 2 )f n 
m N 2 m L n - q l 
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Substituting in Eq. (4.4) we obtain the simplified Goldberger-Treiman relation 

G mA (q 2 )f„ = 2m N Cj(q 2 ) (4.6) 

in complete analogy to the well known GT relation which holds in the nucleon sector. Pion pole 
dominance therefore fixes completely the ratio C^(q 2 )/C^(q 2 ) as a pure monopole term 



Ci(q 



in 



N 



m 2 ~q 2 ■ (47) 

The goal here is to calculate C^(q 2 ), C^{Q 2 ) and G n N&(Q 2 ) and check the GT relations using 
dynamical DWF. The relevant three-point functions required for the calculation of these FFs are 
obtained at a minimal extra cost using the sequential propagators produced from the optimized 
nucleon to A source S\ and in addition S2 which is also used for the electromagnetic transition 
study of the subdominant FFs. The detailed expressions are given in Ref. 
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Figure 5: The ratio C£/C^ versus Q 2 . The dotted 
line refers to the DWF results and is the pion pole 
dominance prediction of Eq. (4.5). The solid line is a 
fit to a monopole form. 



Figure 6: Q -dependence of the pseudoscalar tran- 
sition form factor G^na- The solid line is a fit to 
pion pole dominance form of Eq. (4.9). The dashed 
line is a linear fit. The strong coupling constant 
gxNA is the value at Q 2 = 0. 



In Fig. ^| we compare our results for Cg using DWF to those obtained previously using the hy- 
brid action and quenched Wilson fermions at similar pion masses []6[ g]. The Q 2 dependence is well 
described by a dipole Ansatz yielding Cg (0) = 0.970(30) and a dipole mass m A = 1.588(67) GeV. 
This is to be compared with the value m A = 1.28 ±0.10 GeV extracted by a dipole fit to the avail- 
able experimental data [10]. As in the case of Gm\(Q 2 ), we observe a flatter slope for the lattice 
data, reflected in the larger value of the axial mass m& extracted for the lattice results. 

In Fig. |5| we show the ratio Cg /Cg . The dotted line shows the pion pole dominance prediction 
of Eq. ( |P| ) where for m^ and m n we use the lattice values calculated for DWF. The predicted 
curve does not describe the data at low Q 2 i.e. in the regime where strong pion cloud effects are 
expected. Fitting to a monopole form cq/(Q 2 jm 2 + 1) describes satisfactorily the ratio yielding a 
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heavier mass parameter m than the lattice value of the pion mass. This behavior has been observed 
also for the other actions [g]. 

The pseudoscalar form factor G^a(<? 2 ) is determined optimally from the source Si with a 
pseudoscalar current operator insertion: 



5f(q;75)= EnC(0,-q;r 4 ;75) 



(7=1 



E N + niN 



fam 2 



6m N 2m q (m 2 +Q 2 ) 



G KNA (Q 2 ) (4.8) 



We use the value f % = 0.1052(7) GeV for the pseudoscalar pion decay constant determined in 
Ref. [|9p. The quark mass m q is calculated through the Axial Ward Identity by constructing a suit- 
able ratio of local-smeared and smeared-smeared two-point functions of the axial and pseudoscalar 
currents [g]. This requires only knowledge of the axial current renormalization Za, which is deter- 
mined to be Za = 0.7197(9) (Yamazaki et al in [jl]]), where also Zy = Za holds up to a small 0(a 2 ) 
error for a chiral action [^]. 

In Fig. H we compare results on G^a(<? 2 ) using dynamical DWF to those obtained with the 
hybrid action and in the quenched theory Dm. The solid line is a one-parameter fit to the form 



Q 2 /ml+\ 



(4.9) 



(e 2 K+i) 2 (e 2 /™ 2 +i) 

expected if the validity of Eq. ( |4.7| ) is assumed. The fit-parameter K provides an estimate of the 

strong coupling g^NA = Gwa(0) = 9.6(2). A straight line fit of the form G k n A {Q 2 ) ~ ( 1 —A^j 

as shown by the dashed line, would lead to an estimate g n NA = 13.9(6). Thus a reliable evaluation 
of g X NA requires further understanding of the behavior at low Q 2 and in particular of the decrease 
observed in the hybrid action at Q 2 close to zero. 
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Figure 7: The ratio f K G KNA {Q 2 ) /m N C^(Q 2 ) as a Figure 8: The ratio m N f n G KNA {Q 2 ) /2(m 2 c + 
function of Q 2 relating to the GT validity. Q 2 )C^(Q 2 ) that relates to the validity of Eq. (4.5). 

In Fig. |7| we show the ratio f K G n ^ A {Q 2 ) /m^C^{Q 2 ), which should be unity if the non-diagonal 
GT relation of Eq. ( |4~6| ) is satisfied. Deviations from this relation are evident in the low Q 2 regime 
and they are present for all actions to the same degree which is surprising since one might have 
expected a better behaviour for DWF. At higher momentum transfers (Q 2 > 0.5 GeV 2 ) the relation 
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is satisfied for all actions. On the other hand, the relation given in Eq. ( ft~7| ) that assumes pion pole 
dominance to relate Cg to Cj is satisfied excellently by the lattice data for all three actions. This 
agreement is shown in Fig. || where the ratio mNf K G K NA(Q 2 )/2(m 2 + Q 2 )C^(Q 2 ) is everywhere 
consistent with unity. 

5. Summary and Conclusions 

The nucleon to A electromagnetic, axial and pseudoscalar transition form factors are calculated 
using Nf = 2 + 1 dynamical domain wall fermions for pion mass of 0.33 GeV. The dominant 
form factors Gm\ and Cj show slower falloff with Q 2 as compared to experiment. A possible 
explanation maybe that the pion cloud is still not fully developed, at pion mass of 0.33 GeV. We 
examined the Goldberger-Treiman relations and found that they are satisfied for Q 2 > 0.5 GeV 2 as 
was previously observed for Wilson fermions and when using a hybrid action. Pion pole dominance 
relating the axial form factor Cg and the pseudoscalar form factor G k na is satisfied for all values of 
Q 2 irrespective of the lattice action used. Extraction of the strong coupling constant gxNA requires 
special care since we need a better understanding of the low Q 2 behavior of the pseudoscalar matrix 
element. A calculation on a finer lattice using domain wall fermions is underway to check for any 
cut-off effects as well as obtain results on the subdominant and phenomenologically interesting 
electromagnetic quadrupole form factors. 
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